NMR and dc-susceptibility studies of NaVGe206 
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We report the results of measurements of the dc magnetic susceptibihty x a-nd of the ^^Na nuclear 
magnetic resonance (NMR) response of NaVGe2 06, a material in which the V ions form a network of 
interacting one- dimensional spin S = 1 chains. The experiments were made at temperatures between 
2.5 and 300 K. The xiT) data suggest that the formation of the expected low-temperature Haldane 
phase is intercepted by an antiferromagnetic phase transition at Tat = 18 K. The transition is also 
reflected in the '^'^Na NMR spectra and the corresponding spin-lattice relaxation rate T^^{T). In the 
ordered phase, T^^{T) decreases by orders of magnitude with decreasing temperature, indicating 
the formation of a gap of the order of 12 K in the magnetic excitation spectrum. 

PACS numbers: 75.10.Pq, 75.50.Ee, 76.60.-k 
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1. Introduction 

Quasi-one-dimensional (ID) spin systems have recently 
been the subject of both theoretical and experimental 
investigations. R was shown theoretically [1] that for 
chains of antiferromagnetically coupled integer spins the 
ground state is separated by an energy gap Ajj from the 
magnetic excitations. This was later confirmed experi- 
mentally for many systems [2], including the inorganic 
compounds AgVPzOg [3], CsNiClg [4], YBaNiOg [5], as 
well as for organic substances [6] . 

In 1999, a new member, LiVGe206, was added to the 
list of quasi-lD magnets [7]. From the temperature de- 
pendence of the magnetic susceptibility, the expected 
spin gap A/z/Zcb was estimated to be about 19 K. In- 
stead of developing the Haldane phase at low tempera- 
tures, however, the compound exhibits a transition to a 
three-dimensionally (3D) magnetically ordered state at 
25 K. This was established by NMR [8] [9], as well as by 
neutron-diffraction experiments [10]. 

Since many aspects of the magnetically ordered state 
of LiVGe206 remain a mystery, and the V-V interaction 
depends on details of the V-O-V bonds, we consider it as 
worthwile to carefully investigate the magnetic properties 
of other materials related to LiVGe206. In this work, we 
report a study of the related compound NaVGc206. Al- 
though LiVGe206 and NaVGe206 exhibit many similar 
features, they are not identical at all. In particular, the 
differences in the ionic radii of Li and Na seem to be re- 
flected in subtle structural changes. Our dc-susceptibility 
measurements indicate a magnetic phase transition at 18 
K. We provide a detailed analysis of the magnetic suscep- 
tibility of NaVGe206 (and also of LiVGe206), including 
computer-based model calculations. This leads to esti- 
mates of the intra- and interchain coupling parameters, 
indicating that for NaVGe206, the intrachain coupling 
is weaker and the ID-character less pronounced. From 
the estimated values we conclude that the non negligi- 
ble interchain coupling is the main cause for the absence 
of the expected formation of the Haldane state at low 



temperatures. The results of subsequent ^^Na NMR- 
experiments, including the mapping of NMR-spectra and 
measurements of both the spin-lattice and spin-spin re- 
laxation rates, confirm the existence of a phase transition, 
most likely to a 3D antiferromagnetically ordered state. 

The article is organized as follows. In Section 2 we pro- 
vide information on the sample. In Section 3 we report 
on the dc-susceptibility mesurements, and in Section 4, 5 
and 6 we describe and discuss the results of the NMR ex- 
periments. In Section 7 we compare the magnetic proper- 
ties of the related compounds NaVGe206 and LiVGe206. 



2. Crystal Structure and Sample 

NaVGe206 crystallizes with a monoclinic structure, 
space group P2i/c. The NaVGe206 powder sample was 
prepared following the procedure described in [11]. A sto- 
ichiometric mixture of Ge02 (Aldrich, 99.998%), V2O3 
(Aldrich, 99.99%), and Na2Ge205 was pressed into pel- 
lets and annealed at 900° C for 70 h in an evacuated silica 
tube. Na2Ge205 was synthesized by heating a mixture 
of Na2C03 (Aldrich, 99.995%) and GeOz at 760°C for 15 
h in air. According to X-ray powder diffraction data, the 
sample was of single phase with the pyroxene structure. 
The corresponding lattice parameters are a=9. 9600(5), 
6=8.8444(4), c=5.4858(2) A and /3 = 106.50(2)°. The 
crystal structure is shown schematically in Figure 1. It 
consists of chains of isolated, slightly distorted VOg oc- 
tahedra joined at the edges. These chains are linked but 
also kept apart by double chains of distorted Ge04 tetra- 
hedra. 

Taking into account the most likely oxidation states 
of 0^~, Na^+ and Ge^+, the V ions are expected to be 
trivalent. The magnetic moments are thus due to the 
two 3(i electrons localized at the ions, which form 
an S* = 1 system. Assuming that the interaction beween 
these moments is due to exchange interaction mediated 
by the 2p-electrons of the O ions, the intrachain coupling 
J is anticipated to be much larger than the interchain 
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FIG. 2: Magnetic susceptibility x as a function of tempera- 
ture T, measured in a magnetic field of = 2 T. The solid 
line represents the best fit to the data for temperatures above 
100 K, assuming x{T) = Xo + C(T - e)-\ 



FIG. 1: Schematic representation of the crystal structure of netic moment (in units of /ig) is thus 



NaVGe206. The crystalline axes are labelled. The V ions 
are located in the centers of the dark-grey VOe octahedra, 
while the Ge''+ ions occupy the centers of the light-grey Ge04 
tetrahedra. The Na+ ions are represented by dark spheres. 



coupling J± because the former involves only one O site, 
whereas the latter is spread over two of them. 

For our measurements we used a powdered sample, 
with randomly oriented grains and a mass of 127 mg, 
corresponding to 4.04x10^'' mol. The diameters of the 
powder grains were less than 0.1 mm. 



3. Magnetic Susceptibility 

In order to confirm the above mentioned expected oxi- 
dation state of the V ions, we measured the dc magnetic 
susceptibility xiT) between 4 and 340 K, in various fixed 
magnetic fields between 0.01 and 2 T, using a commercial 
SQUID magnetometer. For temperatures above 100 K, 
the data can reasonably be fitted by a Curie- Weiss type 
of law, i.e., x(T) =xo + C{T - 6)"^- The constant term 
Xo is attributed to magnetic background. In Figure 2 
we display an example of measured in a magnetic 

field of = 2 T. In this field, xo = 1-7 x 10^^ emu/mol. 
The paramagnetic Curie temperature is = — 39 ± 5 K, 
which signals the tendency of the V moments to couple 
antiferromagnetically. The value of C is 



Moff 



> 3kB 



C = 2.67 ±0.05 



(2) 



per V ion. This value is close to the one expected for the 
V"^"*" (and not V'*^) configuration, with quenched orbital 
moments of the two 3d electrons [12]. 

The maximum of xC^) is reached at the temperature 
Tmax = 25 K. Below this temperature, a kink in x(T') is 
observed at Tn = 18 K, as demonstrated in the inset of 
Figure 3. The anomaly is much more evident by plotting 
dxl dT{T). The transition temperature T^r is only little, 
if at all, affected by the value of the magnetic field up to 
5 T. 

As already anticipated in Section 2, the chemical com- 
position and the crystal structure of NaVGe206 suggest 
that the V-moments form a quasi-lD 5=1 spin sys- 
tem. Nevertheless, the temperature dependence of the 
susceptibility does not follow the expectations for a col- 
lection of independent S = 1 spin chains. This is evident 
if one compares the experimental data with the results of 
quantum Monte Carlo simulations [13], as shown in Fig- 
ure 4. The simulation was performed for a single chain 
considering the Hamiltonian 



i/ = J^5.-5 



i+l I 



(3) 



C = 0.89 ± 0.04 emu • moh 



K 



(1) 



where the indicated error range takes into account the 
values found for other magnetic fields. The effective mag- 



where Si denotes a spin-1 operator at the i-th site, and 
J is the intrachain coupling. J was adjusted such that 
the maximum in x{T) occurs at 25 K, as is experimen- 
tally observed. At high temperatures the experimental 
and calculated data coincide. Near 25 K, however, the 
experimental susceptibility turns out to be smaller than 
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FIG. 3: dx/dT as a function of the temperature T for dif- 
ferent applied fields. The inset shows x(^) at H = 2 T. 
The dotted line indicates the low-field transition temperature 
Tn = 18 K. 



expected for the isolated S = 1-chains (see Figure 4) , sug- 
gesting the existence of some non negligible interchain 
coupling. The simulated and experimental x(T) data 
also differ significantly at low temperatures. In partic- 
ular, the continuous decrease to zero of the calculated 
revealing the influence of the Haldane gap, is not 
observed experimentally. Instead a kink in x(^) at 18 K 
indicates a magnetic phase transition. 

As is well known from series expansions [16], even 
a very small interchain coupling Jj^ {J±/J > 0.026) 
leads to the quenching of the Haldane gap and induces 
thrcc-dimonsonal antiferromagnetic order (see also [17]). 
As indicated by x(T), this indeed seems to occur in 
NaVGe206, suggesting a non- negligible interchain cou- 
pling Jj_ in this compound. 

From the measured values of Tmax and Tn we can es- 
timate the magnitudes of both the intra- and interchain 
coupling, and thus their ratio, as follows: Assuming that 
Tmax is unaffected by a small interchain coupling, quan- 
tum Monte Carlo simulations of an isolated chain provide 
a relation between J and T^ax of the form 



J 



1.32 



±0.02 



(4) 



The value of the interchain coupling can be estimated 
from the transition temperature Tn, using a random- 
phase approximation (RPA). The relevant equation is 



J±/J 



1 



zx'[kBTN/J] 



(5) 



where denotes the staggered spin susceptibility of an 
isolated 5=1 spin chain. Within conventional RPA [18] 




+ S=1 chain, quantum Monte Carlo simulation 



100 200 
7(K) 



o 0.010 
E 











(b) 




















mmnm 










********,^ 




















o ^ 








r™, = i.32j/kB 




o c7 












Tn 


= 0.95 J /kg 













10 20 30 40 

r(K) 



FIG. 4: (a) Experimentally measured xiT) ~xo as a function 
of temperature in i/ = 2 T, and x{T) resulting from quantum 
Monte Carlo calculations for an S = 1 spin chain. Panel (b) 
emphasizes the range T < 50 K. 
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FIG. 5: Interchain coupling as a function of the critical tem- 
perature, according to equation (5). The extrapolation to 
ksTN / J = gives J±/J = 0.02. The paramotors J±/J and 
ksTN / J for NaVGe206 arc indicated by the dashed lines, 
while the dotted lines represent J±/J = 0.1 and the corre- 
sponding UbTn/J- The inset shows the staggered spin sus- 
ceptibility x" of a S = 1 spin chain as a function of tempera- 
ture, as calculated by quantum Monte Carlo simulations. 



z = A equals the number of nearest neighbour chains. 
Recently, it was found that for both the S = 1/2 case and 
the classical Hciscnbcrg model (corresponding to S* ^ 
oo), an improved estimate of J^_/ J results by using a 
renormalized value of 

z = 2.78 (6) 
for J_l/J < 0.2 [19]. We apphed this modified RPA 
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also to the S = 1 case. The result for {J±/J) versus 
[ksTt^/J] is shown in Figure 5, along with an inset for 
the staggered susceptibility x^iksT/J] of an S* = 1 spin 
chain, as obtained from quantum Monte Carlo calcula- 
tions. The limit [/csT/v/J] — > determines the critical 
value J±/J = 0.02, i.e., the interchain couphng neces- 
sary to suppress the Haldane gap, which is in reason- 
able agreement with the value estimated in [16]. For 
J_i/J = 0.1 wc find fcsTjv/J = 0.65 ± 0.02 (see dotted 
lines in Figure 5), a value which is in good agreement with 
the result ksTi^ / J = 0.65 ± 0.01 of an explicit quantum 
Monte Carlo calculation for interacting chains. 

The dashed lines in Figure 5 represent the parameters 
J_l/J and ksT/J for NaVGe206. Using equations (4) 
and (5) for T^ax 25 ± 0.5 K and Ta, = 18 ± 0.5 K, we 
obtain for NaVGe206 

J//CB = 18.9 ±0.5 K (7) 

and 

J^/J = 0.18 ±0.01 , (8) 

i.e., 

Jj_/fcB = 3.4±0.2K . (9) 
4. NMR Spectra 
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FIG. 6: ^^Na NMR spectra of the central line in the paramag- 
netic phase at three different temperatures m H = 7.0495 T. 
The integrated echo intensity (multiplied by the temperature 
T) is represented as a function of the frequency /. The solid 
lines represent the intensity /(/), obtained using Equation 
(20) and setting Ao = 1140 G. The dashed curve (shown only 
for T = 80 K) corresponds to the calculation with Ao = G. 
The dotted line marks the reference position of the ^'^Na line 
(79.389 MHz). 



With our NMR experiments we probed the ^"^Na nuclei 
with a nuclear spin / = |, a gyromagnetic factor of 7 = 
7.0746 X 10^s~^T~-'^ and a quadrupolar moment \e\Q = 
0.108|e| X 10^^** cm^, where e is the electron charge. 

^•^Na NMR spectra were obtained by monitoring the 
integrated spin-echo intensity as a function of the irra- 
diation frequency, in a fixed magnetic field H or, alter- 
natively, as a function of the magnetic field, at a fixed 
irradiation frequency. The spin echo was generated with 
a two-pulse 7r/2-delay-7r spin-echo sequence, irradiating 
a frequency window of about 20 kHz. 

In Figure 6, we show examples of the central line of 
the spectra, corresponding to the 5 — 5 transition, 
measured at three different temperatures above T^r, in a 
fixed magnetic field oi H = 7.0495 T. The background 
echo intensity, due to quadrupolar wings, i.e., the transi- 
tions ±1 ±i, is subtracted. 

The paramagnetic character of the high-temperature 
phase is reflected in Figures 7 and 8. In Figure 7, we plot 
the relative shift with respect to the unshifted frequency 
fo, 



peak 



/peak fo 

To 



(10) 



of the frequency /peak with maximum echo intensity at 
different temperatures, as a function of the Curie- Weiss 
part of the dc magnetic susceptibility x{T), measured at 
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FIG. 7: Relative frequency shift A'peak of the maximal in- 
tensity NMR line as a function of the dc-susceptibility x, for 
temperatures above Tjv. The solid line is the best linear fit 
to the data at ii" = 7.0495 T, above 60 K. 



the corresponding temperatures and in a magnetic field 
oi H = 2 T. The data for the two different fields collapse 
onto the same curve. Above Tmax, K(x) can reasonably 
well be fitted with a linear function, indicating that the 
temperature dependence of the shift of the NMR-line at 
maximum intensity is of purely magnetic origin. In Fig- 
ure 7, we display the best fit to the data at H — 7.0495 
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FIG. 8: FWHM of the echo intensity of the c:ontraI line, 
divided by the ^^Na reference frequency /o, as a function of 
the dc-susceptibility x, for temperatures above Tjv. The solid 
and the dashed lines are the best linear fit to the data at = 
7.0495 T and H = 4.5533 T respectively, for temperatures 
above 60 K. 
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FIG. 9: The solid line represents the shape P"^^^ of the central 
line calculated assuming dipolar coupling between V moments 

and Na nuclei, as a function of the hyporfinc coupling A. The 
peak value ^p'J^j^ and the estimated width aj'' arc indicated. 
The dashed line is obtained by broadcuiug the solid line with 
AA = 132 G (corresponding to the same brodaning Af in 
frequency space as used for T = 80 K in Figure 6) . 



T and for T > 60 K with the function 



-^peak — ^^peakX ~l~ '^peak 



(11) 



Within experimental uncertainty, Cpeak ~ 0. The hyper- 
fine coupUng corresponding to the peak signal is [21] 



^peak = apeak ' Nij,b = 1300 ± 50 G 



In Figure 8, we plot 



ax = 



FWHM 



fo 



(12) 



(13) 



where FWHM is the linewidth and /o is the unshifted 
frequency, at different temperatures, as a function of the 
Curie- Weiss part of the dc-magnetic susceptibility x(T), 
which was measured at the corresponding temperatures 
and in a magnetic field of = 2 T. For T > Tmax a 
linear relation aK{x)j i-^-, 



(14) 



is observed. In Figure 8, the solid line represents the best 
fit to the data for H = 7.0495 T and for T > 60 K. The 
corresponding anisotropic hyperfine coupling is [21] 



aA = tta ■ Nfli 



560 ± 50 G 



(15) 



Prom the value of for H = 7.0495 T we find a fre- 
quency broadening A f = Cala = 14 ± 2 kHz, which is 
attributed to the nonvanishing irradiation width that in- 
deed is of the same order of magnitude. The dashed line 
is the best fit to the data for H = 4.5533 T, keeping Oo- 
fixed, resulting in Af = c^/o = 18 ± 2 kHz, again of the 
same order of magnitude. 

The linear relation between ctk and x breaks down 
for T < Tmax, indicating that below that temperature 



a substantial degree of antiferromagnetic correlations in- 
terferes. 

To complete the analysis of the features of the central 
line in the paramagnetic phase, we compared the exper- 
imentally determined lineshape with the results of sim- 
ulations. They were performed by first calculating the 
lineshape P^^p{A) expected for a purely dipolar coupling 
between the V moments and the Na nuclei, i.e., according 
to the formula 

P'^'P{A)= [ d<j) [ de smO 5{A- A'^''P{e,(l))) . 
Jo Jo 

(16) 

Here, 6 and (j) are the spherical angles describing the ori- 
entation of the V moments of 1 /Lis with respect to the 
cristalline c-axis, A'^^P{9,(f)) is the component of the in- 
duced hyperfine-field at a Na site parallel to the V mo- 
ments and S is the delta-function. Equation (16) im- 
plicitely assumes that (a) the powder sample consists of 
randomly oriented grains; (h) the V moments align along 
the applied external magnetic field. 

The integral in Equation (16) was calculated by ap- 
proximating it with the sum over a 200x200-points- 
lattice g C [G,7r] x [0,27r], 

P^'^A)^ sin0rAA(A-A<i'P(0,</.)) , (17) 

{0,4>)eg 

where Faa is a Gaussian function of width chosen as 
AA w 20 G. The resuhing hneshape P'^'p{A) is repre- 
sented as a function of the hyperfine coupling A in Figure 
9. The peak value is reached at 



4?ak = 125±10G 
while the mean square deviation is 
al'P = 520 G . 



(18) 



(19) 
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While crj^ is in good agreement with the experimental 
value (TA, ^pcak ^ factor of ten smaller than Apeak (see 
Eq.(12)), indicating an additional, relevant non-dipolar 
isotropic hyperfine-field coupling. 

For a fixed temperature T and the corresponding sus- 
ceptibility Xi signal intensity /(/) as a function of 
the irradiation frequency / was obtained by 



/(/)=/0,T- [P^'P(A(/)-Ao)]^^. 

where A and / are related by 

/-/o A 



X 



(20) 



(21) 



Jo,T is a (temperature dependent) prefactor, and [. . .]a/ 
indicates the broadening in frequency space with a Gaus- 
sian function of width A/. 

In Figure 6 we plot the calculated lineshapcs with 
A/ = 15±1 kHz, Ao = 1140 G, and Io,t such that the ex- 
perimental intensity at the frequency /peak is matched. 
The parameter Aq takes into account the above men- 
tioned non-dipolar coupling, whose value has been ad- 
justed to reproduce the experimental data, i.e.. 



Apeak « ^pcak + "^0 



(22) 



The experimental lineshape is clearly very well accounted 
for by the simulations. The deformation range of the 
central NMR-line due to second order quadrupolar effects 
is given by [21] 



144/c 



/(/ + l)-4 



(77^ + 22r; + 25) , (23) 



where fc is the Larmor frequency of the ^■^Na central 
line, I'Q = e^qQ/2h is the quadrupolar frequency of ^^Na 
in the electric field gradient eq and i] is the asymme- 
try parameter of the electric field gradient (eq — Vzz is 
the largest eigenvalue of the electric field gradient ten- 
sor, and 77 = {Vxx — Vyy)/Vzz)- As we shall see below, 
i^Q = 1.25 MHz and t] = 0.4. With / - 79.6 MHz this 

results in Afg' ~ 14 kHz. The effect is more evident at 
high temperatures, because it is less masked by magnetic 
broadening. The small deviations of the order of 10 kHz, 
observed in particular at 150 K, are attributed to such 
quadrupolar effects. 

As shown in Figure 10, with decreasing temperatures 
but still above T/v, we identify features near the cen- 
tral line which are temperature and field independent. 
They are the signals of the wings due to the transitions 
±1 ^ ±i. This claim is supported by the results of 
comparisons of the measured signals with the signal cal- 
culated assuming that the intrinsic part of the central 
line, at frequency fc, is unaffected by the orientation of 
the grains. The shift of the ±| ^ ±i-transitions is given 
by [22] 



[3 cos'^ 9 — 1 — T] sin^ 9 cos 2(f 



(24) 
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FIG. 10: ^^Na NMR spectra around the central line in 
NaVGe206 at two different temperatures above Tjv, and in 
two different fields. The solid lines represent the calculated 
lineshapes due to quadrupolar effects. 



Here, 9 and are the spherical angles describing the ori- 
entation of the principal axes of the electric field gradient 
tensor with respect to the direction of the applied mag- 
netic field. The solid lines in Figure 10 represent the 
calculated signal without the contribution of the central 
line. The best coincidence is achieved with a quadrupolar 
frequency is vq = 1.25 MHz and an asymmetry param- 
eter T] — 0.4. The central frequencies fc were chosen 
to coincide with the values of the average frequencies of 
the central line at the corresponding temperature, i.e., to 
51.38 MHz (at 4.5533 T) and 79.57 MHz (at 7.0495 T), 
respectively. The largest eigenvalue of the electric field 
gradient tensor at the Na-sites is 



Vzz = eq = 



2hvQ 
eQ 



9.58- 10^"V/m 



(25) 



Next, we focus our attention on the ^^Na NMR spectra 
around and below the critical temperature Tn = 18 K. In 
Figure 11, we show examples of such spectra, measured 
in the temperature range between 12 K and 20 K, in 
fixed magnetic fields H = 4.5533 T and = 7.0495 T, 
respectively. 

The phase transition at T/v = 18 K is reflected by a 
drastic change of the lineshape. We note that the central- 
line intensity decreases rapidly between 18.5 K and 16.5 
K, indicating that the paramagnetic phase disappears 
upon crossing the critical temperature from above. At 
the same time, a much broader signal develops below 
18.5 K. At 14.5 K the signal has a well defined, "nearly 
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FIG. 11: ^•■^Na NMR spectra of NaVGc206. The echo inten- 
sity is represented as a function of the frequency / for selected 
temperatures below 20 K. The measurements were performed 
in fixed magnetic fields H = 4.5533 T (upper panel) and 
H = 7.0495 T (lower panel). 



rectangular" shape, which remains unchanged down to at 
least 4 K, as evidenced by the signal obtained from a field 
sweep at a fixed frequency of 80.8 MHz, which is shown in 
Figure 12. The observed width and the evolution of the 
lineshape below Tjv are both similar at different values of 
the magnetic field, suggesting that the broadening of the 
line is due to ordered magnetic moments. This allows, in 
particular, to exclude the formation of the Haldane phase 
below Tjv. Indeed, our NMR results cannot be reconciled 
with the expectations for a non-magnetic ground state, 
where a sharp narrowing of the NMR line and a shift to- 
wards the reference frequency /o should be observed at 
temperatures T < Ah/Zcb ~ 0.4 J/ks = 7.6 K [23]. 

The deviation of the observed low-tcmpcrature line- 
shape from a rectangle, which is expected for a purely 
magnetic broadening due to a fixed internal field, in par- 
ticular the feature in the range between the marks Li and 
Ri in Figure 12, is attributed to first order quadrupolar 
effects. It is not a trivial task to predict the lineshape 
due to the combined effects of an external field, an in- 
ternal field and an electric field gradient. The difficulty 
lies in the uncertainty of the directions of the ordered 
moments and of the principal axes of the the electric 
field gradient tensor at the Na-sites, respectively. There- 
fore, we discuss this aspect only qualitatively. We assume 
that at low-temperature (T ^ T/v) the signal due to the 
5 — 5 transition extends from line L2 to line R2 in 
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FIG. 12: Example of a ^^Na NMR spectrum of NaVGe206 
at 4 K. The echo intensity is plotted as a function of the 
magnetic field H. The measurements were performed at a 
fixed irradiation frequency of 80.8 MHz. 
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Ferromagnetic 
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Ferromagnetic 
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Antiferromagnetic 
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c-axis 


Ferromagnetic 


1720 




Antiferromagnctic 


1720 



TABLE I: Internal field Ji/int at the Na-sites, for different 
orientations of the V-moments and for either ferromagnetic 
or antiferromagnetic interchain ordering. Ordered moments 
of 1.89 //B were assumed. 



Figure 12. The intensities in the ranges between L3 and 
L2, as well as between R2 and R3, reflect the quadrupolar 
wings. Their extension, corresponding to about 1 MHz, 
is compatible with the previously calculated quadrupolar 
frequency vq. With these assumptions, the magnitude of 
the internal field, indicated in Figure 12, is ilint = 1470 
G. 

As a plausibility check for the value of the internal field 
ilint established above, we computed the dipolar field at 
the Na-sites due to thrcc-dimensionally ordered magnetic 
moments at the V sites. The (/-factor for the V'^^ ions 
is given by the effective moment found in equation (2) 
as g = Mcff/v/'5'(S'+ 1) = 1.89. In Table I, we show the 
values of the calculated dipolar fields which are identical 
for all the Na-sites. They were calculated assuming or- 
dered moments of 1.89 /is, oriented along different axes, 
and for either ferromagnetic or antiferromagnetic inter- 
chain order. The obtained values are of the same order 
of magnitude as the experimental result, but they do not 
allow to extract rigorous information about the orienta- 
tion of the ordered moments, except that antiferromag- 
netic order along the a-axis seems very unlikely. The 
uncertainity is, in particular, due to the fact that in the 
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FIG. 13: ^^Na spin-lattice relaxation rate Tf^ as a function 
of temperature T. The maximum of T^^ at the critical tem- 
perature is indicated by the dotted vertical line. The inset 
shows an example of the echo intensity as a function of the 
recovery time; the solid line is the best fit using the recovery 
function (26). 
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FIG. 14: ^•■'Na spin-lattice relaxation rate Tj as a function 
of inverse temperature . The solid line represents the 
function Tf ^ = Cexp(-A/T) with A = 12.5 K. 



relation time r, i.e., wr -C 1. In this scenario [22] 



calculations we considered only dipolar coupling. In ad- 
dition, the moments need not be aligned along the same 
direction, and their magnitude might also be reduced (as 
was suggested in the case of LiVGe206, see [9]). 



5. NMR spin-lattice relaxation rate 

In order to probe the low-energy spin excitations in 
NaVGe206 at low temperatures, we mesured the spin- 
lattice relaxation rate Tj~^(T) by monitoring the recov- 
ery of the ^"^Na nuclear magnetization after the appli- 
cation of a long comb of rf pulses, in the temperature 
range between 2 K and 300 K. The experiments were 
performed by irradiating a frequency window of about 
100 kHz. Above Tjv, the irradiation frequency was cho- 
sen to be in the center of the central line, while below Tm 
we centered the window at approximatively the center of 
the broadened line. In both double recovery 



(27) 



m{t) 



1 - D{0.6e + o.4e -^i 



(26) 



characteristic of / = 3/2 nuclei [24], with t the time vari- 
able and a constant £) » 1, was observed. 

In Figure 13, we display the temperature dependence 
of the spin-lattice relaxation rate. A prominent peak in 
Tj~^(T) reflects the phase transition at = 18 K. Above 
Ttv, T^^ varies only weakly with temperature and is not 
much affected by changing the field from 4.5533 T to 
7.0495 T. Those features are typical for a spin-lattice re- 
laxation driven by the spin-flips of the V'^^ paramagnetic 
moments whose dynamics is characterized by a short cor- 



where 7 is the gyromagnetic ratio of the Na- nuclei, A is 
the hypcrfine coupling and S is the spin of the magnetic 
moments. With Tf ^ = 50 7 = 7.0746 x lO^s-^T"! 
and A 0.1 T, we estimate the correlation time of the 
V^+-moments to be r 1.25 x 10~^^s. 

Below T/v, Tj~^ decreases by orders of magnitude with 
decreasing temperature (a factor of 10^ from 18 K to 
2.5 K, as shown in Figure 14), indicating the formation 
of a gap in the spectrum of magnon excitations in the 
antiferromagnetically ordered state. This is supported 
by T^^{T) well below T/v, which varies according to 

Tfi(T)~Ae-T , (28) 

as is emphasized by the solid line in Figure 14. The 

magnon processes which lead to spin-lattice relaxation 
may be very complicated by involving several real and/or 
virtual magnons. At low temperatures, however, it is 
fair to assume that a two-magnon process, namely, the 
scattering of a magnon at the nucleus with a nuclear 
spin-flip (Raman process) dominates [25]. In this case, 
Tj~^(T) is given by Equation (28), and A corresponds to 
the gap in the magnon spectrum. Such a gap is usually 
attributed to an easy-axis single-ion magnetic anisotropy 
term in the Hamiltonian, of the form D'^^ Si^zSi^z with 
D <0 [26]. 

In our case, the uncertainty of the data at the lowest 

temperatures allows only for an estimate of the value 
of the gap A/ks, of the order of 12 K. If a small 
amount of magnetic impurities, not indicated by the low- 
temperature susceptibility, were present, the actual value 
of the gap could be somewhat larger. 
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FIG. 15: ^•■^Na spin-spin relaxation rate T2 as a function 
of temperature T. The inset shows an example of the echo 
intensity as a function of the decay time; the solid line is the 
best fit using the exponential decay function in equation (29) . 



6. NMR spin-spin relaxation rate 

Our measurements of the spin-spin relaxation rate Tj"^ 
as a function of temperature were made in the range be- 
tween 2 K and 300 K. The spin-echo hfctimc T2 was ob- 
tained by monitoring the spin-echo intensity as a function 
of the pulse delay r in the 7r/2-delay-7r echo sequence. 
The irradiation conditions were chosen to be the same as 
in the Ti-experiments. At all temperatures, the intensity 
decayed according to 



i{t)^ioe n 



(29) 



Here, i = 2{2/it ■ + d) + t^, ^ 2d [28], where 
and are the durations of the pulses, and d is the delay 
between them. The spin-spin lattice relaxation rate was 
then calculated from 



^2 



(30) 
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The last approximation is justified because T2 ^ 

In Figure 15, we display the temperature dependence 

of the spin-spin relaxation rate, measured for H = 7.0495 
T. A prominent peak in T2"^(T) is observed at 17.5 ±0.5 
K, and reflects the phase transition at T^r. Below that 
peak, T2^{T) exhibits an unexpected upturn with de- 
creasing temperature and approaches a constant value at 
temperatures less than 5 K. Generally, is determined 
by the slow fluctuations of the internal magnetic field at 
the Na-sites. Since in our case Tj"^ » TT^, these fluctu- 
ations must basically be given by spin-spin interactions 
between the Na-nuclei. In the antifcrromagnctically or- 
dered state, besides the direct dipolar interaction, also a 





LiVGc2 0(i 


NaVGe2 06 


T^ax (K) 


62 


25 


J/fes (K) 


47.0 


18.9 


Tn (K) 


25 


18 


J±/fcB (K) 


3.4 


3.4 




0.07 


0.18 



TABLE II: Comparison of data extracted from dc- 
susceptibility measurements on LiVGe206 and NaVGe206. 



magnon mediated interaction is expected to play a role. 
At any rate, it is not easy to see how this would lead to 
an increase of T2~^ at lower temperatures and we have no 
reasonable explanation for this particular feature. 



7. Comparison with LiVGe2 06 

From the structural point of view, LiVGe206 is very 
similar to NaVGe206. More precisely, the two substances 
share the same space group P2i/c, and the lattice param- 
eters differ by less than 1.5%. Moreover, both compounds 
contain chains of VOe-octahedra, kept apart by Ge04- 
tetrahedra [7] . Previously, we and other authors reported 
dc-susceptibility and NMR studies on LiVGe206 [7] [8] [9] . 
The results turned out to be qualitatively the same as 
those for NaVGe206 reported in this work. Below we 
focus our attention on the most evident quantitative dif- 
ferences. 

The dc-susceptibility of LiVGe206 X^) shows a max- 
imum at Tmax = 62 K and a kink at T^v = 25 K, from 
which the intra- and interchain couplings may be ob- 
tained in the manner described in Section 3. The cal- 
culated values are shown in Table II. While the inter- 
chain coupling turns out to be approximatively the same, 
the intrachain coupling in LiVGe206 is about a factor 3 
larger than that in NaVGc206. This is somewhat sur- 
prising and indicates that even small differences in the 
orientations of the orbitals involved in the exhange in- 
teraction can affect the effective coupling between the 
moments at the V sites. 

LiVGe206 has been the subject not only of NMR- 
experimcnts, but also of neutron scattering [10] and muon 
spin resonance experiments [29]. It has been concluded 
that the low-temperature phase is antiferromagnetically 
ordered, with a ferromagnetic interchain coupling; the 
magnitude of the ordered moments was reported to be 
1.14/Us, a much smaller value than expected from the g- 
factor inferred from the susceptibility. Using the above 
value for the ordered moment and single crystal NMR- 
experiments results, the authors of [9] suggested that the 
ordering occurs along the c-axis. Our data are consistent 
with a similar scenario for NaVGe206. 
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8. Conclusions 

NaVGe206 can be considered as a quasi l-dimensional 
5=1 spin system only at temperatures much higher 
than the corresponding equivalent of the intrachain cou- 
pling J/fcs ~ 19 K. The expected Haldanc phase, with 
a Haldane gap An/ks of the order of 8 K, does not 
develop because the system orders antiferromagnetically 
below 18 K. The interchain coupling, J±/kB ~ 3.4 K, 
estimated from the x{T) data, is by far sufficient to ex- 
plain the suppression of the Haldane gap. The decay of 
the spin-lattice relaxation rate T^^{T) ^ exp(— A/fc^T) 
observed at T Tjv, indicates the opening of a gap 
A/ks < 12 K in the magnon excitation spectrum of the 
three dimensional antiferromagnetically ordered system. 
This phenomenon is attributed to fluctuations of the two 
3d-electrons orbitals of the ions, leading to an easy- 
axis single-ion anisotropy term in the Hamiltonian. 

Further experiments with single crystals, aiming at ob- 
taining more information on the orientation and the mag- 
nitude of the components of the electric field tensor, the 
internal field and the eventual magnetic anisotropy seem 
to be in order. 
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